Nuclear matter in the chiral quark soliton model with vector mesons 
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We study the nuclear matter solution in the chiral quark soliton model coupled to p and lo vector 
mesons based on the Wigner-Seitz approximation. It is shown that the vector mesons stabilize 
the soliton at high-density region. As a result, the saturation property and incompressibility are 
significantly improved. 
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I. INTRODUCTION 

The idea of investigating dense nuclear matter in 
the topological soliton models has been developed over 
decades. It was first applied for the nuclear matter sys- 
tem vifith the skyrmion centered cubic (CC) crystal by 
Klebanov 0. This configuration was studied further by 
Wiist, Brown and Jackson to estimate the baryon den- 
sity and discuss the phase_transition between nuclear 
matter and quark matter y. Goldhabor and Manton 
found a new configuration, body-centered cubic (BCC) 
of half-skyrmions in a higher density regime -3]. The 
face centered cubic (FCC) and BCC lattice were stud- 
ied by Castillejo et al. Q and the phase transitions be- 
tween those configurations were investigated by Kugler 
and Shtrikman Recently, the idea of using crystal- 
lized skyrmions to study nuclear matter was revived by 
Park, Min, Rho and Vento with the introduction of the 
Atiyah-Manton multi-soliton ansatz in a unit cell 

The soliton model incorporating quark degrees of 
freedom into each soliton was also considered in 
80's. Achtzehnter, Scheid and Wilets investigated the 
Friedberg-Lee soliton bag model with a simple cubic lat- 
tice ■ Due to the periodicity of the background poten- 
tial, the solution of the Dirac equation has the form of 
the Bloch waves, "ipkif) — e''^'^0fc(r) where </)fc satisfies 
the same periodic boundary condition as the background 
potential. 

The Wigner-Seitz approximation was used for the anal- 
ysis of the crystal soliton model with quarks. In this 
ansatz, a single soliton is placed on the center of a spher- 
ical unit cell. Then the lowest energy level ("bottom" of 
the band) for the valence quarks becomes s-state. The 
appropriate boundary conditions at the cell boundary 
should be imposed on the quark wave functions as well 
as the chiral fields. This simple treatment sheds light on 
the nucleon structure in nuclear medium. Soliton matter 
within this approximation have been extensively studied 
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by using various nucleon models such as the the chiral 
quark- meson type model ^ | ^ ITol ITlj . Friedberg-Lee 
soliton bag model dl. l3.ll3lll4j |. the Skyrme model .1^ . 
The non-zero dispersion of the lowest band [ll| and the 
quark-meson coupling [l^ were also examined within this 
approximation. 

The chiral quark soliton model (CQSM) can be inter- 
preted as the soliton bag model including not only va- 
lence quarks but also the vacuum sea quark polarization 
effects explicitly [H [13 El [l^- The model provides 
correct observables of a nucleon such as mass, electro- 
magnetic value, spin carried by quarks, parton distribu- 
tions and octet, decuplet SU{i) baryon spectra "2^, '2l'|. 
Amore and De Pace studied nuclear matter in the CQSM 
using the Wigner-Seitz approximation and observed the 
nuclear saturation |22|. They examined the soliton solu- 
tions with three different boundary conditions imposed 
on the quark wave function. However the obtained sat- 
uration density was lower than the experimental value. 
They thus concluded that such discrepancy is originated 
in the approximate treatment |2,3| | of the sea quark con- 
tribution. 

In Ref . ■ we studied the nuclear matter in CQSM 
and observed splitting of the nucleon-A spectra. The 
vacuum polarization was treated exactly and a relatively 
shallow saturation was obtained. However for the value 
of the constituent quark mass M reproducing the octet 
and decuplet baryon spectra, the soliton breaks even at 
low densities. 

In this paper, we construct the matter soliton solutions 
including tt, ct, p and uj. The role of p, uj is to produce the 
short range effects of the nuclear force and stabilize the 
solution at high densities. It is straightforward to include 
the p meson in the CQSM, but the uj meson requires some 
technique since the Hamiltonian is no longer real. To 
overcome this difficulty, we apply two different methods 
proposed for the free nucleon system H^H^ and compare 
the obtained results. 

This paper is organized as follows. In the following 
section, we present the basic formulation of CQSM with 
vector mesons. Two distinct formulations for solving 
the non-Hermitian eigenvalue problem are reviewed in 
Sec lIIII In Sec lIVI we show how various cutoff param- 
eters and coupling constants are determined within the 
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chiral perturbation regime. In Sec0 the extension of 
the model to the nuclear matter within the Wigner-Seitz 
approximation is presented. The numerical results are 
shown in Sec I VII Sec lVIII is devoted to summary and 
conclusions. 



II. THE CHIRAL QUARK SOLITON MODEL 
WITH p,uj MESONS 

The CQSM was originally derived from the instanton 
liquid model of the QCD vacuum and incorporates the 
non-perturbative feature of the low-energy QCD, spon- 
taneous chiral symmetry breaking (SCSB). The semi- 
bosonized version of the Nambu-Jona-Lasinio model also 
inspires the CQSM model with the SCSB. In these de- 
scription, the Euclidean vacuum functional with vector 
mesons can be defined as EE Ha 



X exp 



(1) 



where — X]a=o ^■''"/-^ i'^^ — 1) ^''^ vector gauge 
fields for the vector mesons. The SU(2) matrix 



{775 



1 + 75 



U 



1-75 



(2) 



with 



U = exp [it ■ (p/fn) = —{a + it ■ -k) . 



describes the chiral fields, a and tt represent scalar 
sigma meson and pseudoscalar pion fields respectively. 
■0 denotes quark fields and M is the dynamical quark 
mass. /,r is the pion decay constant and experimentally 
/tt ~ 93MeV. Since our concern is the tree-level pions 
and one-loop quarks according to the Hartree mean field 
approach, the kinetic term of the pion fields which gives 
a contribution to higher loops can be neglected. Due 
to the interaction between the valence quarks and the 
Dirac sea, soliton solutions appear as bound states of 
quarks in the background of self-consistent mean chiral 
field. Nc valence quarks fill the each bound state to form 
a baryon. The baryon number is thus identified with the 
number of bound states filled by the valence quarks [23| ■ 
The B = 1 soliton solution with only chiral fields has 
been studied in detail at classical and quantum level in 
Refs. Ellllllllllllillll. 

Integrating over the quark fields in Eq.(^, we can ob- 
tain the effective action ScS for the mesons 



^cff ~ Sp + Sm, 

Sf = -iNJndet f - Af 



(4) 
(5) 

(6) 



where Sm is the term derived from the semibosonized 
version of NJL action [2^. Sp can be divided into real 
and imaginary parts: 



Sf — Sr + Si, 

Sr ^^Trlogil^^m, 

sj = ^TTiog{m-^m 



(7) 
(8) 

(9) 



where il^ = y — MU^'^ is the modified Dirac oper- 
ator. After performing the Wick rotation for the Dirac 
operator, i.e. xq = —1x4 and Vq = —iVi, we can obtain 
the one-quark Hamiltonian H with vector fields 

i(3I^ = -dr-H, (10) 
H = a-P + iVi + a-V + (3A4U^^ . (11) 

Here r denotes the Euclidean time. Note that as r, V4 is 
supposed to be Hermitian in Euclidean space, H is now 
non-Hermitian. 

To obtain the B — 1 soliton solution, let us impose the 
hedgehog ansatz on the chiral field 



U{r) = exp{iF{r)f ■ t) 

= cos F (r) + if ■ T sin F (r) 



(12) 



The only possible ansatz for the isoscalar-vector field lj^ 
realizing zero grandspin is the one whose spatial compo- 
nents vanish (uji = 0), 



(13) 



(3) Parity invariance requires the isoscalar-axialvector me- 
son field V4 to vanish in the static limit. Note that 
corresponds to the physical lo meson. For the isovec- 
tor and vector meson fields let us impose the spherically 
symmetric ansatz 

1 



(14) 



where the indices a,i and fc run from 1 to 3. corre- 
sponds to the physical p meson. The boundary conditions 
of F(r) for the B = 1 soliton solution are given by 



F{0) = -TT, F{oo) = 0. 



(15) 



Regularity requires the following boundary conditions for 
Lu and G, 



w'(0) = 0, tj(oo) = 0, 
G(0) = 0, G{oo) = 0. 



(16) 



Substituting these ansatz into Ea. ljlll) . one obtains the 
effective Hamiltonian 

H = a ■ p + iijj{r) + - (a x f) • TG{r) 

+l3M{cosF{r) + ij5T • rsini^(r)). (17) 
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As stated above, H is non-Hermitian since the real func- 
tion of Lo{r) makes H complex- valued. The eigenvalue 
problem of H is solved using the method developed by 
Kahana-Ripka ([l^l, see also Sec lV|l . 

Once the eigenvalue of H , e^, is obtained, the eigen- 
values A„ ^ of the operator dr + H (fTn|l are determined 
by 



= —iQr 



(18) 



where irj„ = i{2n + l)7r/T with (n = 0,±1,±2,..) and 
e^, are the real and imaginary part of eigenvalues of 
the Hamiltonian H17|l . The quark determinant is ex- 
pressed in terms of the eigenvalues A„.p as 



Sr = \11 log(A„,^A;^), Si = \Y. log 



A. 



71, /i 



.(19) 



Since the real part Sr diverges as log for large mo- 
menta p, we apply the proper time regularization [26j |. 
The real part of the sea quark energy from can be de- 



rived from e"^^ ~ e ^s 



E. 



El 




(20) 



Similarly the imaginary part of the sea quark energy is 
derived from e'^^ ~ e~'^™c^ as 



(21) 



where 




(22) 



The static energy for the vector mesons is given by 



where gi and 52 a-re coupling constants determined in 
the subsequent section. The total energy i?tot is defined 
by the sum of these energies plus (three times of) valence 
quark energy (see Sec lIII|l . Field equations for the meson 
fields can be obtained by demanding that the total en- 
ergy be stationary with respect to variation of the profile 
function, 



= 



(24) 



where (p denotes any of the meson profile F,G ov uj, which 
produces 

Sir) sin F{r) = P{r) cos F{r) , 



—-fopix,x), 

ihoTor ■ t))p{x,x) 
— ((7 X r) • t)p{x,x) 
w[r) = g2Nctr / —b(x,x) 



(25) 



where p{x,y), b{x,y) are the quark scalar density and 
the quark number density respectively. Their specific 
forms are presented in the next section. 



III. THE EIGENVALUE PROBLEM FOR THE 
NON-HERMITIAN HAMILTONIAN 

In order to solve the eigenvalue problem with the non- 
Hermitian Hamiltonian, let us first introduce the left and 
right eigenstate 



I.e. 



(26) 



with the normalization condition (^'^ | ^'^) ~ (5^^. For 
convenience we separate the Hamiltonian into Hermitian 
and non-Hermitian part 



H ^ Hi 



e 



lOJ 



(27) 



where Hq is the Hermitian part. The fact that Hq and 
ljj are both Hermitian implies | \f =1^*^. There are 
two distinct ways 0, to extract the physical spectra 
from the eigenequations H26|l . 

In Ref . |25l| . the Wick rotation from Euclidean to 
Minkowski space has been performed to the Hamiltonian. 
Since the time component of the vector fields becomes 
iujQ the eigen equations are reduced to 



Era^\^J d^'xGirf + \^J dPxLoirf . (23) Defining 



one can write the valence quark energy as 



(28) 



(29) 



(30) 



where rii, is the valence quark occupation number. Then 
the total energy is given by 

Etat [F, UJ, G] — i?val + -E^TOc + ^vac + ^m- (31) 

By substituting l|31|l into Ea. (|24(l . we obtain the equation 
of motion (|25|l where the quark scalar density and the 
quark number density are given by 



p{x,y) — F^-ipip'^+F -ip ijj 
b{x,y) = F+^il;+ - F-^-%p- 



(32) 
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FIG. 1: The total energy as a function of in the formalisms 
of Refs.|25ll and |2^ . In the formalism of Ref.|25||. the solution 
does not exist for g^j > 3.8. guj ~ 4.6 is predicted in the chiral 
perturbation analysis (see Sec.5). 




Unfortunately, this simple method do not produce so- 
lutions for the predicted value of the uj meson coupling 
constant g^^ required in the chiral perturbation analysis. 
In Fig. n we show the total energy of the soliton as a 
function of . As can be seen, the soliton survives only 
up to gui ~ 3.8 whereas the chiral perturbation analysis 
predicts the value around g^^ ~ 4.6. 

In the second method, Eq. H26(l is solved directly. Then 
the real and imaginary part of the one particle energy 
eigenvalue are derived as [26l |. 

= 1{{^;\H\^,) + {^,\H\^;)) 

^ {^^ I He - {^i I He 

+ {^l I He \^^) + I He (34) 

where we employed the decomposition 1^'^) =1^*^) + i \ 
vp^) and (^f* 1= (VP^ I +i(\E'^ |. The valence quark 
energy is given by the same form as in the first method 

E^i = N,^n,\e% E^, = iV, ^ n,sign(e^)|6^ |. (35) 

If V 



FIG. 2: The binding energy in the formalisms of Refs.|2!Tl| 
and [2^ . where g^—2.7 is used for the former and g^ ~ 4.6 
for the latter. 

The total energy in Euclidean space is 

i?tot[F, u;, G] = + + + eU + Er.. (36) 

Thereby the total energy in Minkowski space is inter- 
preted as 

Etot [F, U), G] = £'^1 + -E4i + -E'rac + 'E'vac + E,n- (37) 

The equation of motion takes the form in H25|) with the 
replacement of p{x, y) and p{x, y) by 

p{x,y) = pl^' + pr + pT + Pr, (38) 

pl^\x,y) + pT{x,y) 
= [^^{x)^^{y) - i^iixmy)] (n. + f^), 

pT\x,y) + pr{x,y) 
= [i'^^ixWAy) + i^i{x)^!}{y)] {n, + fl), 

(39) 

b{x,y) = b]^' + bY' + bY + bT" , (40) 
bt{x,y) + bT{x,y) 

= Y [i^^{xW{y) - ^i{xWJ{y)\ (n. + fl), 

bT\x,y) + bnx,y) 

= -J2 [^!^ix)^iHy) + ^^(=r)^f (y)] {n, + /«), 

1/ 

(41) 
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FIG. 3: The profile functions of tlie cliiral fields for 
i?=0.75,l,1.5,2fm and the free {R oo) solution. 



where fn and // are the imaginary part of the regurahzed 
action written by 



= sign(e^)AA,- 
= sign(6^^)AA,. 



1 



(ei/A)cxp(-(6^/A)2), 



IV. THE COUPLING CONSTANTS OF THE 
VECTOR MESONS 



In this section, we derive the coupUng constants for 
the vector mesons. Expanding the real part of the ef- 
fective action (jHl up to second order in the heat-kernel 
expansion ^2^], one obtains 



.fl,(2) 



167r2 



d^xT 0, 



A" j 



ctr Ul'd^Ud^U^ + \f^^,FI:^ + \Fl,Fr\ , (42) 



where 



We renormarize the vector meson fields V., 



-iguiV^ for later convenience. Then the 



total Lagrangian is given by 



tr 



where 



0^ ^2 



F^. = d^v^ - d^v;: + 5p[v;^ 



(43) 



(44) 



Comparing Eq. (|43|l with the massive Yang-Mills La- 
grangian, the following relations for the parameters are 
obtained, 



f 



9p = 



A/2 



6AP 2 6Af2 

9^ = 



2 ' 



4ffi' " 4ff2 



4/^ 

2 



(45) 
(46) 
(47) 



The experimental values are /^=93MeV, M=350MeV, 
Mp=770MeV, M^=783MeV, and hence 

A = 649MeV, gi = 1.39, 32 = 0.34, 

gp = 9.22, g^ = 4.61. (48) 

As a reference, let us note that in Ref. p9|. 

gp = 8.32 , g^ = 4.41 (49) 

are adopted for the MIT bag model. 



V. THE NUCLEAR MATTER FORMULATION 

In this section we describe the numerical method of 
eigenproblem of the Hamiltonian (|17l) . The Hamiltonian 
with hedgehog ansatz commutes with the parity and the 
grandspin operator given by 

K = j + T/2 = l + a/2 + r/2, 

where j,l are respectively total angular momentum and 
orbital angular momentum. Accordingly, the angular ba- 
sis can be written as 

m)KM) = Cmr. (E ClX^^rn)\U))\U) . 



(50) 
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FIG. 4: The profile functions of p meson for i?=0.75,l,1.5,2fm 
and tiie free {R oo) solution. 



FIG. 5: The profile functions of uj meson for 

i?=0.75,l,1.5,2fm and the free {R oo) solution. 



For B = 1 solution, following states are possible 



\0) = \{K K + -)KM), 
\l) = \iK K-^)KM), 
\2) = \{K+IK+^)KM), 
\3) = \iK-lK-^)KM). 

With this angular basis, the normalized eigenstates of the 
free Hamiltonian in a spherical box with radius R can be 
constructed as follows: 



u 



V 



V 



'■KM 

(b) 
KM 

(a) 

KM ■ 

ib) _ 
KM 



--Nk 
--Nk 



^ejK+i{kr)\2) 

i^tjK{kr)\l) 
-^7jK-iikr)\3) 

'i^tjK+i{kr)\2) 
-cj~jK{kr)\0) 

iuj+jK-i{kr)\3) 
^7JK{kr)\l) 



(51) 



with 



Nk = 



^R^'^jK+iikR) 



-1/2 



(52) 



•^i>0'^e,<o = sgn(efe),w,^>o,w+<o = k/{ek + M). 
The u and v correspond to the "natural" and "unnatural" 
components of the basis which stand for parity (—1)^ 
and (—1)^+^ respectively. The momenta are discretized 
by the boundary conditions jK{kiR) = 0. The orthogo- 



nality of the basis is then satisfied by 
drr^jKihr)jK{kjr) 



R 



drr'^jK±iihr)jK±iikjr) 



sAjK±i{k^R)r 



(53) 



Let us examine the boundary conditions for the Dirac 
and chiral fields to construct the nuclear matter solu- 
tion in the Wigner-Seitz approximation. When the back- 
ground chiral fields are periodic with lattice vector a, the 
quark fields would be replaced by Bloch wave functions 
as ip{r + a) — e^''''^ifj{r). In the Wigner-Seitz approxi- 
mation, however, the soliton is put on the center of the 
spherical unit cell with the radius R {a = 2R) and the 
dispersion k is assumed to be zero. Then, R is related to 
the baryon density through the relation 



3 



PB 



47ri?3' 



(54) 



For the Dirac eigenstates, modification in the basis is 
needed. For odd number of K , the boundary condition 
is same as the free case with 

jKihR)^0. (55) 

For even K, the following conditions must be satisfied 



jK+iik^''^R) ^0, for u^km^'^k'm^ 
jK-iikf'R) = 0, for u^^j,!, v''j^\i . 



(a) Ja) 



(56) 



From the conditions (|56|l together with the equations of 
motion (|25|) . we find the following boundary conditions 
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FIG. 6: The "upper" u{r) and the "lower" w{r) positive 
component of valence quark wave functions for various cell 
radius R with the boundary condition w{R) = 0. 



FIG. 7: The negative component of valence quark wave func- 
tions. The label is as same as Fig. |S| 



for the profile function F{r) 

ct'(O) = (7'{R) = 
7r(0) = tt{R) = 



F{Q) = -TT,F{R)=Q, (57) 



whieh guarantees the periodicity and the unit topological 
charge inside the cell. Also, for vector meson profiles 
w(r), G(r) we find the conditions instead of Eg. ljltill 



tj'(O) = 0, uj'iR) = 0, 
G{0) = 0, G{R) = 0. 



(58) 
(59) 



We solve (|25ll selfconsistently with boundary conditions 
(|56|l - (|59|l for varying R, from infinity (isolate the soliton) 
to origin (infinite density matter). 



VI. THE NUMERICAL RESULTS 

Since the equations of motion for mesons H25|l and the 
Dirac equations for quarks H26|l are highly non-linear, we 
solve these equations numerically by selfconsistent anal- 
ysis. Using trial profiles for tt, p and uj mesons which 
satisfy the appropriate boundary conditions, we solve 
the Dirac equation (|26|l . From Ea. H25|l . the profile func- 
tions F{r), G{r) and uj{r) are uniquely determined by 
using the eigenstates of Ea. H26|l . The new profiles pro- 
duce new eigenstates. These procedures are repeated un- 
til self consistency is attained. We chose he quark mass 
M = 350MeV which is mostly used and turned out to be 
the best choice for obtaining the various experimental ob- 
servables [23, |^ ■ We performed the computation from 
i?=5 (infinity) upto the value where the soliton solution 
breaks. The most dense solution is obtained for i?=0.5fm 
which corresponds to the density p=1.91fm~'^ ~ Hpjv- 



Note that /9Ar=0.17fm is the standard saturation den- 
sity. 

In Fig. 121 the binding energies computed in the two 
methods are shown. In the first method, the binding 
energy is too large, and the solutions disappear at rela- 
tively lower densities. Therefore, we focus our attention 
to the numerical results obtained in the second formalism 
hereafter. 

Fig.Oshows the self-consistent profile functions for free 
{R ^ oo) and for the various values of the cell radius R. 
Figs. 21 |S1 show meson profile functions of G{r),uj{r) re- 
spectively. Fig. O shows the real part of the quark wave 
functions. Nonvanishing values of the upper component 
at the cell boundary u(i?) come from the zero-mode ele- 
ments in the basis. The imaginary part of the quark wave 
functions derived in the second formulation is shown in 
Fig.EI 

In Fig. we present the results of the total energy of 
the soliton as a function of R, for various values of the 
meson couplings. For only the pion coupling, the satu- 
ration can not be observed and the soliton disappears at 
very low density. Adding the p meson, the soliton sur- 
vives at higher density, but no saturation is observed. In 
the case of tt — tj, the total energy is enhanced and the 
soliton disappears at low density. The behavior is simi- 
lar to the case of the pion coupling only. The saturation 
can be observed only when all the meson couplings are 
incorporated. The saturation point is at i?=1.0fm cor- 
responding to the density p — 0.24fm~^ and the binding 
energy 53MeV. The density is very close to the experi- 
mental value for nuclear matter which is i?=l.lfm corre- 
sponding to the density p = 0.17fm~^. The binding en- 
ergy is deeper compared to its empirical value, 16MeV. 
In Fig. 1^1 we show the binding energy for varying the 
constituent quark mass, which is the only free parameter 
in our model. The soliton survives for the constituent 
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FIG. 8: The total energy of soliton for the various combina- 
tions of the mesons. 



FIG. 9: The total energy varying the constituent quark mass 
M. 



quark mass in the range SOOMeV < M < 440MeV. As 
M increases, the saturation point moves to the lower 
density and the binding energy becomes shallower. For 
M=350MeV, the solution reaches to the highest density. 

Let us estimate the nuclear incompressibility since it 
gives an important information for the saturation prop- 
erty of the matter. In Ref. the authors studied the 
soliton matter in the Friedberg-Lee model with quark- 
meson coupling using the Wigner-Seitz approximation. 
They estimated the incompressibility K with the formula 



K = R 



2 (PEb 



(60) 



and obtained K ~ llTOMeV. The experiment predicts 
ii:=100-500MeV and generally K ~ 200MeV. In our 
previous analysis l24j with only tt mesons, we obtained 
K ~ 400MeV. Our new result predicts K ~ 270MeV. 



VII. SUMMARY 

In this paper we studied the nuclear matter solutions in 
the chiral quark soliton model. We adopted Wigner-Seitz 
approximation and investigated the saturation property 
of the matter solutions. To improve the qualitative be- 
havior at the saturation point, we introduced the p, w 
mesons in the model. The p meson can be incorporated 
in a straightforward manner. On the other hand. Incor- 
porating the Lo meson makes the Hamiltonian complex- 
valued and requires some technique. We tested two dis- 
tinct formulations proposed in Rcfs. 25, 26J to solve the 
non-Hermitian eigenvalue problem. It turned out that 
the latter method produces the better behavior of solu- 
tions especially at high-density region. We performed 
the chiral perturbation analysis to determine the me- 
son properties, i.e., the coupling constants and the cut- 



TABLE I: We present the p, uj couphng constant gp, g^, the 
total energy of the free nucleon E^^t (its the experimental 
value is 939MeV), the binding energy Eb (16MeV), the sat- 
uration density ps (0.17fm~^) and the incompressibility K 
(210±30MeV) for the constituent quark mass M=350, 375, 
400MeV. 



M[MeV] 


9p 


g^ 


£;f^r[MeV] 


£B[MeV] 




K[MeV] 


350 


9.22 


4.61 


1095 


53 


0.24 


271 


375 


9.88 


4.94 


1070 


36 


0.13 


242 


400 


10.54 


5.27 


1044 


29 


0.11 


274 



off parameter. We found that stable nuclear matter 
solutions exist when tt, p, uj mesons are included with 
330 <M< 420 [Me V]. 

From Fig. |S1 one can speculate that the attractive 
property of the p meson makes the total energy smaller 
and allows the solution to survive at the high density 
regime. The repulsive property of the lo meson makes 
the total energy larger and creates a short range (high 
density) core. Thus, the model having both effects can 
produce a stable nuclear matter. 

Although our formulation realizes a improved nuclear 
saturation property, the binding energy is still large. 
Throughout the calculation, we set the constituent quark 
mass M=350MeV, in which the solution survives at high- 
est density. By varying the value of M in a few ten MeV, 
the saturation property should be improved (see Fig. 1^. 
TABLE m shows the computed saturation properties for 
Af=350, 375, 400MeV. M=375MeV seems to be the best 
choice, but the binding energy is still large. Another at- 
tempt to improve the binding energy may be to introduce 
heavier mesons, like axial vector meson ai (tt, a, p, uj and 
ai are the mesons which do not vanish for hedgehogs 
|25j|'). It is known that ai meson has attractive property 
like p meson and then it would make the matter softer. 
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Taking account the Fermi motion would also make the We arc grateful to Kouichi Saito for fruitful discussions 
saturation energy shallower. and comments. 
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